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We present recent results of the rare semileptonic decay B —> XsBl~ ■ We particnlarly focus on 
higher order electroweak corrections in the Standard Model (SM) as well as 0[as) corrections 
in Two-Higgs-donblet models (THDM), both of which are computed within an effective field 
theory approach. The calculation of higher order electroweak corrections reveals the presence 
of enhanced electromagnetic logarithms \n{mf,/m\) in the differential branching ratio. The 
inclusion of C>(as) in the THDM reduces the scale dependence of the corresponding Wilson 
coefficients significantly. 


1 Introduction 

The inclusive decay B with I = e or pt probes the Standard Model directly at the 

one loop level and is therefore sensitive to new physics. Its branching ratio has been recently 
measured by both Belled and BaBar^^. The experimental results for the differential branching 
ratio, dB{B —> XJ~^l~) / dq‘^, integrated over the low dilepton invariant mass region, 1 GeV^ < 
< 6 GeV^, read 

BiB^Xsl^r) = (1.493 ± 0.504l°;^^J) X 10“® (Belle), (1) 

BiB^XJ+r) = (1.8 ± 0.7 ± 0.5) X 10"® (BaBar) , (2) 

leading to a world average 

B{B^Xsl+l-) = (1.60 ± 0.51) X 10"® . (3) 

Other appealing features of the decay B —> Xsl~^l~ are on the one hand the possibility to 
obtain complementary information compared to the less rare decay B Xgj, on the other 


hand, precision data on both the experimental and theoretical side can be achieved. Indeed, 
the experimental errors in the branching ratio are expected to be substantially reduced in the 
near future. On the theoretical side, the predictions are quite well under control because the 
inclusive hadronic B —> decay rate for low dilepton mass is well approximated by the 

perturbatively calculable l!.~ decay rate. Thanks to the recent (prac¬ 

tically) complete calculationof the Next-to-Next-to-Leading Order (NNLO) QCD 
corrections, the perturbative uncertainties are now below 10%. 


However, at the leading order in QED, the branching ratio is proportional to ae^(^), giv¬ 
ing rise to a ±4% scale uncertainty when the renormalization scale of Oem is changed from 
^ = O(mfe) to /r = 0{Mw)- This uncertainty can be removed by calculating higher order 
electroweak corrections. The auth ors of Ref.calculated the QED corrections to the Wilson 
coefficients. In a recent calculation^^ we confirmed the results of Ref.and, in addition, com¬ 
puted corrections to the differential branching ratio that originate from QED matrix elements 
of four-fermion operators. It turns out that the latter corrections are numerically relevant since 
they are enhanced by large electromagnetic logarithms ln(m^/m|), which originate from these 
parts of the QED bremsstrahlung corrections where the photon is emitted collinearly by one of 
the outgoing leptons. 


In the task for the search for new physics beyond the Standard Model it is not only relevant to 
obtain precise predictions for the process in question with the Standard Model as the underlying 
theory, but also to perform precision calculations for extensions of the Standard Model (SM). 
In many extensions of the SM, there are additional one-loop contributions in which non-SM 
particles propagate in the loop. If the new particles are not considerably heavier than those of 
the SM, the new contributions to these decays can be as large as the SM ones. One should try 
to get information on the parameters in a given extension - here the two-Higgs-doublet models - 
from all processes which allow both a clean theoretical predictio n and an acc urate measurement. 
This means that precision studies similar to those for B —> where higher order 

QCD corrections are crucial, should also be done for the process B —> Xsl^l~. We focus on 
QCD corrections to the Wilson coefficients Cq and Cio in two-Higgs-doublet models (THDM), 
which we have calculated in Ref.^^ and which prove to decrease the scale dependence of the 
Wilson coefficients significantly. Diagrams with neutral Higgs-boson exchange are neglected. 
This omission is justified in the type-H model, if the coupling parameter (m^/M^y) tan/3 is 
sufficiently smaller than unity. In this case the operator basis is the same as in the SM. Only 
the matching calculation for the Wilson coefficients gets changed by adding the contributions 
where the flavor transition is mediated by the exchange of the physical charged Higgs boson. 
While th^ e extra pieces are known for the coefficients Cj, Cg and Cio to two-loop precision 
|12 1 1 j | 1^ 112] £q]. quite some time, the corresponding re sult s for Cg have been first calculated in 
Ref.^^ and were confirmed and first published in Ref.l^. 


2 Theoretical framework 

2.1 Effective Theory 

To describe decays like B —> Xsl^l~ we use the framework of an effective low-energy the- 
opyi mpdl^l l with five quarks, obtained by integrating out the heavy degrees of freedom. In 
the present case these are the t-quark, the and boson as well as — for the THDM 
calculations — the charged Higgs bosons whose masses Mh are assumed to be of the same 
order of magnitude as M^. We only take into account operators up to dimension six and set 
nig = 0. In these approximations the effective Lagrangian relevant for our application (with 
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^ Ci{^i)Oi + Y, CiQ{fi)OiQ + Ck{^i)Ob . 

.2 = 1 2=3 
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The operators Oi{^) in the first sum are needed for both our SM and THDM calculation. They 
read 


AGf 

^eff ^qcdxqed(ui d, s, c, b, e, /r, t) H ^ts^tb 


Ol = {SL1^J'^CL) {cl7^T%l) , O 2 

O 3 = {sLlfjbL)Y.q{Ql^Q) , C>4 

O 5 = {sLli,lulphL)Y.q{qi^YYQ) , C>6 

O7 = ^mb{sL<Tf^'^bR)F^, , Os 

O 9 = $isLlpbL)El{in), OlO 


{SLlpCL){cLl^bL) , 

{sL7pT%L) EqiQl^T'^Q) , 

{sLlplulpT^^bi) , (5) 

j-m,{sLa^'^T%R)Gl, , 
^{sL7pbL)El{h'"l5l), 


where T“ (a = 1,...,8) are the SU{3) colour generators, and gs and e are the strong and 
electromagnetic coupling constants, respectively, q and I appearing in the sums run over the 
light quarks {q = u ,..., b) and the charged leptons, respectively. 

Once QED corrections in the SM are considered, five more operators need to be taken into 
account. They can be chosen as 


O 3 Q = {sLlpbL) J2q Qqiqi^q), 

O4Q = isLlpT%L) E, Qqiqj^T-q), 

05Q= {sL7pi'yp2lp3bL)J2qQq{q7^^'y^^'y^^q), (6) 

Ogq = isL7pl7p2l^l3'd^''bL)T,qQqiq7^^7^^7^^T‘^q), 

Ob= J2 [isL'ypi'yp2lt^3bL){h^^'y^^7^^b) - iisL7pbL)ih^b)] . 

where Qq are the electric charges of the corresponding quarks (| or — g). 

The Wilson coefficients Gi{g,) are found in the matching procedure by requiring that con¬ 
veniently chosen Green’s functions or on-shell matrix elements are equal when calculated in 
the effective theory and in the underlying full theory up to O [(external momenta and light 
masses)^/M^], where M denotes one of the heavy masses like or Mr. The matching scale 
Hw is usually chosen to be at the order of M, because at this scale the matrix elements or 
Green’s functions of the effective operators pick up the same large logarithms as the correspond¬ 
ing quantities in the full theory. Gonsequently, the Wilson coefficients Gi{gw) only pick up 
“small” QCD corrections, which can be calculated in fixed order perturbation theory. 


2.2 Two-Higgs-doublet models 

In the following we consider models with two complex Higgs-doublets 4>i and (/> 2 . After sponta¬ 
neous symmetry breaking these two doublets give rise to two charged (H^) and three neutral 
{H^, h^, Higgs-bosons. When requiring the absence of flavour changing neutral currents at 
the tree-level, as we do in this paper, one obtains two possibilites, the type-I and the type-II 
THDmI 22], xhe part of the Lagrangian relevant for our calculation is the Yukawa interaction 
between the charged physical Higgs bosons and the quarks (in its mass eigenstate basis): 


Note that there are several normalizations of O 7 - C9io on the market. 








The couplings X and Y are 

X = —cot/3, Y = cot/3 (type-I), 

X= tan/3, Y = cot/3 (type-II), 

where tan /3 = U 2 /ui, with and V 2 being the vacuum expectation values of the Higgs doublets 
(/>i and cl) 2 , respectively. 


3 Standard Model corrections to the decay B —> Xsl~^l 


3.1 Electromagnetic corrections to the differential branching ratio 

An important quantity in the studies of the rare decay B —> Xsl~^l~ is the differential branching 
ratio, dB{B —> Xsl~^l~) / dg^, with respect to the invariant mass of the final state lepton pair. 
The differential branching ratio as a function of has a region of on-shell intermediate cc- 
resonances like the J/'k or the 'kh This region is therefore not accessible perturbatively and 
one distinguishes two g^-windows below and above the cc-resonances respectively. The low-g^- 
window is taken to be from 1 GeV^ < < 6 GeV^, whereas the high-g^-window is considered 

for q"^ > 14.4 GeV^. M any properties — advantages and disadvantages — of each window are 
summarized in Ref. We shall restrict ourselves to the low-g^-window here and write the 
differential decay width as 

dT{B Xs l^l-) _ - sf 


ds 
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where we have introduced the notation s = q"^ fm^. The effective Wilson coefficients contain 


all c orrect ions relevant for the calculation up to NNLO in QCDl^EISlZISIHl and up to NLO in 

QEDEnni. 

The last term contains finite gluon bremsstr ahlu^ s corrections^^. Further more , we 


have included the non-perturbative 0{l/m‘^) correctionsl '^'^'^'* l and 0{l/m‘^) corrections^^ 

In order to minimize the uncertainty stemming from and the CKM angles, we nor¬ 

malize the decay width to the measured semileptonic one. Furthermore, to avoid introduction 
of spurious uncertainties due to the perturbative b XcCu phase-space factor, we follow the 
B ^ Xg'y analysis of Ref.where 


C = 


Vub 


Vcb 


r{B Xcei?) 
T{B ^ X^eh) 


(9) 


was used instead. The f acto r C = 0.58±0.01 has been recently determined from a global analysis 
of the semileptonic data^^. Our expression for the B —> Xs£^i~ branching ratio finally reads 


dB{B XJ+j-) 
ds 


= B{B^Xceu) 


exp 


VtlVtb 


Vcb 


4 4>££(s) 


where <k„ = 1 -|- aem, is defined by 


T{B 


V Gpl^b,pole 12 ^ 
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As stated earlier, the branching ratio has at leading order in QED a ±4% scale uncertainty 
due to the renormalization scale dependence of aem- The removal of this uncertainty requires 


















the inclusion of higher order electroweak corrections, namely QED corrections to the Wilson 
coefficients^^. These corrections have been known for quite a while and imply that the RGB’s 
for the coupling s ar e coupled differential equations that have a perturbative expansion in ag 
and ctem- In Refthe results of Reffor all the two-loop anomalous dimension matrices that 
are relevant for the running of the Wilson coefficients from high scales of order 0{^w) down to 
scales of order 0{fj,b) were confirmed. 

In addition, corrections to the differential br anch ing ratio that originate from QED matrix el¬ 
ements of four-fermion operators were computed^^. The loop corrections are not free of infrared 
divergences and must therefore be considered together with the corresponding bremsstrahlung. 
The dilepton invariant mass differential decay width is not an infrared safe object with respect to 
emission of collinear photons. Hence, QED corrections contain an explicit electromagnetic log¬ 
arithm ln{mf/m‘j), which stems from these parts of the QED bremsstrahlung corrections where 
the photon is emitted collinearly by either of the final state leptons. These enhanced terms al¬ 
ways arise as single electromagnetic logarithms accompanied by an electromagnetic coupling Pem 
and therefore do not get resummed. The log-enhanced parts of the QED matrix elements disap¬ 
pear after integration over the whole phase space available but survive and remain numerically 
important when we restrict to the low dilepton invariant mass region, 1 GeV^ < < 6 GeV^, 

that we consider. Their numerical impact on the differential branching ratio integrated over the 
low-g^ window is about -|-5 % for final state electrons and about +2 % for final state muons. 
The numerical results for the branching ratio integrated over the low-g^ region read 


= (1.59 ± 0.11) X 10 
Bee = (1.64 ± 0.11) X 10 


-6 


( 12 ) 

(13) 


However, the large effect for electrons gets reduced in size roughly to that of muons once the 
experimental resolution for collinear photons is taken into account. Precise nur nber s and more 
profound explanations on the advent of the collinear logarithm are given in Ref.^^. 


3.2 The forward backward asymmetry Afb 

Another appealing quantity of the decay B - 
metry Afb defined as 


Xgl+l- is the so-called forward backward asym- 


Afb{(1 ) = 


2 \_ dSKii/dq^ [cos 9i > 0) — dSKn/dq^{cos 9 1 < 0) 


dSKii/dq‘^ {cos 9 1 > 0) -|- dSKu/dq^ {cos 9 1 < 0) ’ 


(14) 


where 9^ is the angle between the positively charged final state lepton and the initial state B- 
meson in the restfra me o f the final state lepton pair. The forward backward asymmetry is also 
nicely reviewed in Ref.l^. As it is defined as a difference of two quantities over the corresponding 
sum, it is almost insensitive to hadronic uncertainties since the la tter tend to cancel in the ratio. 
In the SM, the forward backward asymmetry has a zero at^-' . 


ql = (3.76 ± 0.22theory ± 0.24^J GeV^ 


(15) 


which is also subject to receive contributions from large electromagnetic logarithms!^. 

The branching ratio and the forward backward asymmetry of the decay B —> Xsl~^l~ are 
important for yet another reason. Gontrary to the branching ratio of the rare decay B —> Xg^, 
which is at lowest order proportional to both the branching ratio of H —> Xsl^l~ and 


the forward backward asymmetry are sensitive to the sign of . Ghanging the sign of Cj'' 
results on the one hand in a shift of the branching ratio. This shift is so large that the value of 
the branching rati o ge ts moved out of the experimentally allowed range. Therefore the SM sign 


^e// 


of is favored 


ehi 


On the other hand, a change of the sign of removes the presence 



of a zero in the forward backward asymmetry 




Hence already a rough measurement of the 


branching ratio and the shape of the forward backward asymmetry can yield useful information 


about the sign of . 


The determination of the sign of is crucial since it allows to strongly 


constrain the parameter space of certain SUSY models*^ 


4 QCD corrections to the Wilson coefficients Cg and Ciq in the THDM 

For the following it is convenient to expand the Wilson coefficients Ci{^w) as 

CAmv) = cf (w) + + 0(9,'). (16) 

The analytic formulae for the QCD corrections to t he W ilson coefficients Cg and Cio in the 
THDM at the matching scale are given in Ref. In this section we briefly illustrate 
the impact of these two-loop corrections on We then introduce a rescaled Wilson 

coefficient (see Eq. uni) 

^ . 1 47r 

ClQ,Y^{^J^w) = ^ —7 -rC'io,H(A‘w)- (17) 

asiidw) 

In Fig. n we plot the quantities 

and “I > (1^) 


i.e. two approximations of Cio.H as a function of the charged Higgs boson mass Mh for the MS- 
and for the pole mass scheme of the t-quark mass. As inp ut par ameters we use as{Mz) = 0.119, 
^poie _ Q QgY^ — gQ _4 QgY and = 0.231 xhe upper frame shows these 

quantities at the relatively low matching scale = Mw- As in this case and rnt{^w) 

are numerically almost identical, the one-loop approximations (dotted and dashed lines) are 
close to each other. The inclusion of the two-loop corrections, however, considerably lowers 
the (absolute) size of the coefficient for all values of Mh considered. In the lower frame a 
higher matching scale of = 300 GeV is chosen. As in this case and rnt{^w) differ 

considerably, the renormalization scheme dependence of the one-loop results is rather large. 
When taking into account the two-loop corrections (solid and dash-dotted lines), the scheme 
dependence is drastically reduced. 

Looking at the renormalization group equation (RGE)13for Ciq.h, one finds that Cio.H does 
not run in QGD, i.e. 

= C'io,H(l^vy), (19) 

where the low scale fih is of the order of m^. In Fig. [2 we show the dependence of C'io,H(A*b) on 
the matching scale nw for Mh = 300 GeV. It can be clearly seen that the inclusion of the two- 
loop contributions significantly lowers the dependence on For > 250 GeV, Cw,Y{{^^b) at 
two-loop precision is nearly /rvv-iadependent. For between Mw and 250 GeV the two-loop 
Wilson coefficient varies about ±4%, whereas the corresponding one-loop coefficient varies about 
± 11 %. 


5 Summary 

The rare decay B Xsl^l~ is subject of many contemporary analyses in particle physics since 
it is a promising channel in the search for new physics beyond the SM. We computed NLO 
QED corrections to the matrix elements of effective operators and found that these corrections 
include terms that are enhanced by large electromagnetic logarithms \n{ml/mf) whose numerical 
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Figure 1: Dependence of the rescaled Wilson coefficient C'io,H(/rw) (see Eq. II71 1 on the charged Higgs boson 
mass Mh at the matching scale = Mw (npper frame) and = 300 GeV (lower frame). The dashed (dotted) 
line is the one-loop contribution expressed in MS -scheme (pole-mass scheme) of the t-quark mass, while the solid 
(dash-dotted) line includes the two-loop corrections in the respective scheme. 
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Figure 2: Dependence of the rescaled Wilson coefficient C'io,H(/i6) on the matching scale ^,w (see Eq. 11711 1 for 
Mh = 300 GeV. The dashed line shows the one-loop contribution expressed in MS -scheme for the t-quark mass, 
while the solid line includes the two-loop corrections in the same scheme. 


impact on the low-q^ branching ratio is in the range of several percent. The zero of the forward 
backward asymmetry is also expected to get shifted by this type of corrections. 

We also showed that the inclusion of QCD corrections to the charged Higgs induced con¬ 
tributions to Cio(/tw) in type-I and type-II THDM significantly lowers the dependence of this 
Wilson coefficient on the scale ^w- 
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